clf


zeroone=0

if zeroone
    x0=0
    x1=1
else
    x0=0
    x1=3
end


hold on
c1 = 2*(2-sqrt(2));
c2 = 2;



for g=0:.01:3
    if g<c1
        a1 = (-4 - g + sqrt(16*g+g.^2))./(-4+2*g);
        if t < a1
            u = 
        end
            
    end
end



g = xvals(c1,c2);
a1b = (2*(2 - 2*g + sqrt(-g.^2 + g.^3)))./(4 - 4*g + g.^2);

plot(g, a1b, 'r')

a3b = (-6*g.^2 + 4*sqrt((-1 + g).*(g.^2)) + 2*g.*(2 + sqrt((-1 + g).*(g.^2))))./(((-2 + g).^2).*(-g + sqrt((-1 + g).*(g.^2))));
plot(g, a3b, 'g')

%plot([c1 c1], [0  1]);
%[g' a1a' a3a' a1b' a3b']

g = xvals(x0,x1);
a3c = (-6 - g + sqrt(-28 + 44*g + g.^2))./(4*(-2 + g));
plot(g, a3b, 'k')


%a2c = -(1/2) + a3 + (a3.^2*g)/2

legend('a1a', 'a3a', 'a1b', 'a3b')



